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ABSTRACT 

Smoothing operation to make continuous density field from observed 
point-like distribution of galaxies is crucially important for topological or 
morphological analysis of the large-scale structure, such as, the genus statistics 
or the area statistics (equivalently the level crossing statistics). It has been 
pointed out that the adaptive smoothing filters are more efficient tools to 
resolve cosmic structures than the traditional spatially fixed filters. We 
study weakly nonlinear effects caused by two representative adaptive methods 
often used in smoothed hydrodynamical particle (SPH) simulations. Using 
framework of second-order perturbation theory, we calculate the generalized 
skewness parameters for the adaptive methods in the case of initially power-law 
fluctuations. Then we apply the multidimensional Edgeworth expansion method 
and investigate weakly nonlinear evolution of the genus statistics and the area 
statistics. Isodensity contour surfaces are often parameterized by the volume 
fraction of the regions above a given density threshold. We also discuss this 
parameterization method in perturbative manner. 

Subject headings: cosmology: theory — large-scale structure of the 
universe 

1. Introduction 

The large-scale distribution of galaxies is one of the most important sources to study 
formation and evolution of cosmic structures. Now there are two ongoing large-scale 
redshift surveys of galaxies, the Sloan Digital Sky Survey (SDSS, Gunn & Weinberg 1995) 
and the Anglo- Australian Telescope 2dF Survey (Colless 1998). These surveys will bring 
us enormous information of three-dimensional galaxy distribution and are expected to 
revolutionarily improve our knowledge of the large-scale structure in the universe. 

The two-point correlation function of galaxies, or its Fourier transform, the power 
spectrum are simple as well as powerful tools, and have been widely used to quantify 
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clustering of galaxies (e.g. Totsuji & Kihara 1969, Peebles 1974, Davis & Peebles 1983). 
These two quantities are based on the second-order moments of matter fluctuations. When 
the fluctuations are Random Gaussian distributed, the two-point correlation function or the 
power spectrum contains full statistical information of the fluctuations. Even though the 
initial seed of structure formation is often assumed to be random Gaussian distributed as 
predicted by the standard inflation scenarios (Guth & Pi 1982, Hawking 1982, Starobinsky 
1982), this simple assumption has not been observationally established. Moreover, the 
cosmic structures observed today are more or less affected by nonlinear gravitational 
evolution. Therefore, only with these two quantities, we cannot study the large-scale 
structure properly. 

Other statistical methods have been proposed and is expected to play complimentary 
roles to the traditional analyses based on the second-order moment. The skewness 
parameter characterizes the asymmetry of one point probability distribution function 
(PDF) of the density field and has been investigated in deep (Peebles 1980, Fry 1984, 
Juszkiewicz, Bouchet & Colombi 1993, Bernardeau 1994, Scoccimarro 1998, Seto 1999). 
Beside higher-order moments (such as skewness), there are other statistical approaches 
designed to directly measure the geometrical or morphological aspects of galaxy clustering. 
Connectivity of the isodensity contour is an interesting target for these approaches. For 
example, the genus statistics were proposed by Gott, Melott & Dickinson (1986), the 
area statistics (equivalently the level crossing statistics) by Ryden (1988), and percolation 
analysis by Klypin (1988). In addition, the Minkowski functionals recently attain much 
attention (e.g. Minkowski 1903, Mecke, Buchert & Wagner 1994, Schmalzing & Buchert 
1997, Kerscher et al. 1997). 

To analyze observed cosmic structures, smoothing operation becomes crucially 
important in some cases. The observed galaxies are distributed in point-like manner, but 
geometrical or morphological analyses, such as, the genus statistics, are usually based on 
continuous (smoothed) density field (see also Babul & Starkman 1992, Luo & Vishniac 
1995). We have traditionally used filters with spatially fixed smoothing radius for analyzing 
the large-scale structure. Even though this method is the simplest from theoretical point 
of views, other possibilities are worth investigated. The local statistical fluctuations due 
to the discreteness of particles are determined by the number of particles contained in the 
smoothing kernel. If we use a spatially fixed filter, we can measure smoothed quantities 
at overdensity regions relatively more accurately than at underdense regions. As a result, 
quality of information becomes inhomogeneous. This inhomogeneity is caused by the 
simple choice to spatially fix the smoothing radius. There must be more efficient methods 
to resolve cosmic structure from particle distribution (Hernquist & Katz 1989). Actually, 
Springel et al. (1998) have pointed out that the signal to noise ratio of the genus statistics 
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is considerably improved by using adaptive smoothing methods. Adaptive methods are 
based on Lagrangian description, use nearly same number of "particles" (mass elements) 
to construct smoothed density field (Hernquist & Katz 1989) and are expected to be less 
affected by discreteness of mass elements. Therefore, it seems reasonable that we can 
resolve cosmic structures more efficiently, using these methods. 

In this article, we perturbatively analyzed quantitative effects caused by adaptive 
smoothing methods. We pay special attention to three representative examples, the 
skewness parameter, the genus statistics and the area statistics. As the skewness is basically 
defined by the one point PDF, we can, in principle, discuss it without making continuous 
density field. But its analysis is very instructive to see nonlinear effects accompanied with 
adaptive smoothing methods. As a first step, we mainly study the density field in real space 
and do not discuss the effects of biasing (Kaiser 1984, Bardeen, Bond, Kaiser & Szalay 
1986, Dekel & Lahav 1999). 

This article is organized as follows. In §2 we describe basic properties of the adaptive 
smoothing and introduce its two main approaches. Then perturbative formulas are derived 
for each of them. In §3 we discuss the skewness parameters of the density field smoothed 
by these two adaptive methods. We evaluate them using second-order perturbation theory. 
Some of results in this section can be straightforwardly applied to the density field in the 
redshift space. Then we discuss weakly nonlinear effects of the genus and the area statistics 
using the multidimensional Edgeworth expansion method explored by Matsubara (1994). 
To characterize the isodensity contour, parameterization based on the volume fraction 
above a given density threshold is often adopted. In §4.1 we discuss this parameterization 
in perturbative manner. In §4.2 we explicitly evaluate the generalized skewness which is 
closely related to the genus and the area statistics. In §4.3 and §4.4, we show the weakly 
nonlinear effects on these two statistics with various smoothing methods. We make a brief 
summary in §5. 



The (unsmoothed) density contrast field S(x) at a point x is defined in terms of the 
mean density of the universe p and the local density p{x) as 



2. Adaptive Smoothing Method 



S(x) 




(1) 



P 



In this article we assume that the primordial density fluctuations obey Random Gaussian 
distribution which is completely characterized by the (linear) matter power spectrum. 
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Unless we state explicitly, we limit our analysis in the real space density field. But some of 
our results are straightforwardly applied to the redshift space quantities, as shown in the 
next section. 

Isotropic filters with spatially constant smoothing radius R have been traditionally 
used to obtain continuous smoothed density field 5fr(x) as follows 

5 FR (x) = j dx'5(x')W(\x' - x\;R). (2) 

Here the function W(x; R) is a spatial filter function and the subscript F indicate the fixed 
smoothing. Most of theoretical analyses in the large-scale structure have been based on this 
fixed smoothing method. As for the functional shape of W(x; R), two kinds of functions 
are often used (e.g. Bardeen et al. 1986, Matsubara 1995). One is the Gaussian filter and 
defined as 

W(x; R) = (2n)-V 2 R- 3 exp (-x 2 /2R 2 ) . (3) 
The other one is the top-hat filter and has a compact support as 

W( X ;R) = l 3 J^ (4) 
v ' ; \0 (|au| > R). K J 

In this article we mainly use the Gaussian filter. This filter is useful for quantifying the 
large-scale structure from observed noisy data sets. In addition, algebraic manipulations for 
the Gaussian filter are generally much simpler than for the top-hat filter. 

Next let us discuss the basic properties of adaptive smoothing methods (Hernquist & 
Katz 1989, Thomas & Couchman 1992, Springel et al. 1998). The essence of these methods 
is to change the smoothing radius R as a function of position x according to its local density 
contrast. With a given spherically symmetric kernel W, we determine the smoothing radius 
R(x) so that the total mass included within the kernel becomes constant. 

p J dx'(\ + 5(x'))R(xfW(\x' - x\;R(x)) = pR 3 . (5) 

The radius R(x) becomes smaller than the standard value R in a overdense region and 
becomes larger in a underdense region. In a system constituted by equal mass particles 
as in standard N-body simulations, the smoothing radius R(x) is determined so that the 
total number of particles in a filter becomes constant. Thus adaptive smoothing is basically 
Lagrangian description and their smoothing radii are closely related to the resolution of 
spatial structures. 

We can solve the variable smoothing radius R(x) in equation (5) by perturbatively 
expanding the deviation 8R(x) = R(x) — R. In this procedure we regard the density 
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contrast S as the order parameter of the perturbative expansion. After some calculations 
we obtain the first-order solution as follows 

5R(x) = ~5 FR (x)R + 0{5 2 ). (6) 

This simple result seems quite reasonable with the relation below. 

R(x) 3 (l + 5 FR (x)) = R 3 (l + 0(5 2 )). (7) 

This relation roughly shows that the total mass within the smoothing radius R(x) does not 
depend on position x. 

With the variable smoothing radius R(x) (solution for eq.[5]) we can practice adaptive 
smoothing. As pointed out by Hernquist & Katz (1989) for the smoothed particle 
hydrodynamics (SPH), there exist two different methods (gather and scatter approaches) to 
assign the smoothed density contrast field at each point x. The gather approach is simply 
use the solution R(x) at the point x in interest and the smoothed field is formally written 
as 

5 GR (x) = J dx'S(x')W(\x' -x\;R(x)) -C(R), (8) 

the subscript G indicates the gather approach. In this case, the volume average of the first 
term in the right hand side dose not vanish and we have added a term C(R) so that the 
total volume average of 5gr(x) becomes zero. 

In the scatter approach we use the solution R(x') for each point where a mass element 
exists. We can write down the smoothed field at x as 

5 SR (x) = J dx'5(x')W(\x' - x\;R(x')), (9) 

the subscript S represents the scatter approach. In this case the volume average becomes 
zero. We only dilute the mass element at point x' with the density profile proportional to 
W(\x' — x\; R(x')) around that point. Note that the spatial dependence of the smoothing 
radius R(-) is different between equations (8) and (9). 

Next we evaluate equations for Sgr(x) and 5sr(x) up to second-order of the density 
contrast 5 using perturbative solution of the smoothing radius R(x) = R + 5R(x) given in 
equation (6). The results are given as 

s GR (x) = SM^-^M^R^M^ + lR^l + ois 3 ), (10) 

Ssr(x) = S FR (x) - 1 1 dx'd R W(\x' - x\;R)6(x')5 FR (x') + 0(5 3 ). (11) 
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The formula for the scatter approach is somewhat complicated, compared with the gather 
approach. Also in numerical analysis, the scatter approach requires higher computational 
costs (Springel et al. 1998). This reflects nonlocal character of the smoothing radius. 

Equations (10) and (11) show apparently that the corrections due to the adaptive 
methods start from second-order of 5. Therefore, their effects are expected to be comparable 
to second-order (nonlinear) effects predicted by cosmological gravitational perturbation 
theory (Peebles 1980, Fry 1984, Goroff et al. 1986). Adaptive smoothing methods modify 
the quantities which characterize the nonlinear mode couplings, such as the skewness 
parameter of density field. 

If we use the Gaussian filter, the leading-order correction for the scatter approach is 
expressed as follows 



/ 



dk dl 

(2vr) 3 (2tt) 3 6XP 



(2Z 2 + fc 2 + 2k ■ l)R 2 
2 



(k + 1VR 2 

5(k)8(l) exp[i(fc + I) ■ x] J , (12) 



where 5(k) and 5(1) are the Fourier coefficients of the density contrast and defined as 

5(k) = J dx5(x) exp(—ik ■ x). (13) 

Formula (12) is useful to quantitatively evaluate nonlinear effects caused by the scatter 
approach. 



3. Skewness 

In this section we investigate modifications of the skewness parameter S caused by the 
two adaptive methods. Skewness is a fundamental quantities to characterize asymmetry of 
the one point PDF of the density field (Peebles 1980, Fry 1984, Juszkiewicz, Bouchet & 
Colombi 1993, Bernardeau 1994). It is defined as 

1 /2 

where the angular bracket (•) represents to take the ensemble average and a(= (<5 2 ) v ) is 
the rms fluctuation of 5. Here, we discuss the leading-order contributions for the numerator 
(5 3 ) and denominator a 4 . As we have already commented in §1, the skewness parameter 
can be discussed without making continuous density field. It can be basically defined by the 
count probability distribution function, and spatial relation between one region and another 
one is unnecessary (e.g. Gaztahaga 1992, Bouchet et al. 1993, Kim & Strauss 1998, and 
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references therein, see also Colombi, Szapudi & Szalay 1998). Therefore our effort in this 
article to resolve cosmic structures by using the adaptive smoothing might be irrelevant for 
observational determination of the skewness parameter. But perturbative analysis in this 
section is very useful to grasp nonlinear effects caused by the adaptive smoothing methods 
and become basis for studying statistics of isodensity contours such as the genus statistics 
or the area statistics discussed in the next section. 

The leading-order contribution for the rms fluctuation a is written in terms of the 
linear (primordial) power spectrum (5i(k)Si(l)) = (27r) 3 5D rc (fc + l)P(k) (5\(k) : linear 
mode, Sr)rc(-)'- Dirac's delta function). With the Fourier transformed filter function w{kR) 
we have 

4 = (<£(*)) = / j^y 3 P(k)w(kR) 2 , (15) 

where the suffix R is added to explicitly indicate the smoothing radius R. Throughout in 
this article, we use power-law spectra P(k) as 

P(k) = Ak n , -3<k<l. (16) 

for these scale-free models the normalization factor A becomes irrelevant and we can simply 
put A = 1 below. Furthermore, as shown later, the skewness parameter does not depend 
on the smoothing radius in our leading-order analysis. From equation (15) we have the 
variance a\ for the Gaussian filter as 

2 f°° dk _ k 2 B 2 i2 -n-3 /3 + n\ 

The integral (15) logarithmically diverges for n = —3, but skewness S is well-behaved in the 
limit n — > —3 from above. As it shows interesting behavior at this specific spectral index, 
we also discuss quantities at n = — 3 regarding them as the limit values. 

Calculation of the third-order moment (5 3 ) is more complicated than that of the 
variance a 2 discussed so far. When the initial fluctuation is random Gaussian distributed 
as assumed in this article, the linear contribution for the third-order moment becomes 
exactly zero due to the symmetric distribution of the density contrast 5 around the origin 
5 = 0. Nonlinear mode couplings induce asymmetry in this distribution. Therefore, we 
resort to higher-order perturbation theory. The leading-order contribution for the skewness 
parameter without smoothing operation is given by Peebles (1980) in the case of Einstein 
de-Sitter background as 

S=%. (18) 

It is convenient to use the Fourier space representation to calculate the third-order moment 
for the smoothed density field. Following the standard procedure, we expand a nonlinear 
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Fourier modes of overdensity 8 and the (irrotational) peculiar velocity field V as (Fry 1984, 
Goroff et al. 1986) 

8(x) = S^x) + 8 2 (x) + ■ ■ ■ , 

V(x) = V 1 (x) + V 2 (x) + ---, (19) 

where 8i(x) and V\(x) are the linear modes and 8 2 (x) and V 2 {x) the second-order modes. 
We perturbatively solve the continuity, Euler and Poisson equations, 

^8(x) + -V[V(x){l + 8(x)}\ = 0, 
ot a 

§-V(x) + -[V(x)-V]V(x) + ^V(x) + -V(/ ) (x) = 0, 

V 2 <f){x) -47ia 2 p{t)8(x) = 0, 

where a represents the scale factor. The second-order solution in fc-space is given as 

r dl 

fc(fc) = j T^r-M l M k - l ) J ^ k - z )> ( 2 °) 



or m cc-space 



52{x) = I (w(w ^ kX5im{k " z) J(Z) k " °' (21) 



(2tt) 3 (2tt) 3 
where the kernel J is defined by 



Ak, l) = kl + K) + ^(i + ^)+hl- ^)^#. (22) 



2 V ' 2 U 2 I 2 J 2 V 7 /c 2 / 2 

The factor fT(fi, A) weakly depends on the density parameter Q and cosmological constant 
A as shown in the fitting formula (Matsubara 1995, see also Bouchet et al. 1992) 

K(Q, A) ~ ^Q- 1 ^ - A (1 _ log 10 n) . (23) 
In the ranges of two parameters Q and A 

0.1 < n < 1, 0.1 < A < 1, (24) 

the difference of K(Q, A) from K = 3/7 is within 8%. Therefore, in the following analysis 
we basically study the Einstein de-Sitter background and use K = 3/7. 

Using the second-order solution (21) we can derive the well known formula for the 
third-order moment as follows (Juszkiewicz et al. 1993) 

. r dk dl 

S3r ) = 6 J j2^j2^ mmJ{K l M kR M lR M\ k + l \ R )- (25) 
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Let us simplify this six-dimensional integral dkdl. In the case of the Gaussian filter 

w(kR) = exp(-k 2 R 2 /2), (26) 
we can change (5^) to the following form (Matsubara 1994) 

4) = ^( 5/ 220 + 7/131 + 2/222), (27) 

where we have defined 

POO roo rl 

I abc = dk dl duexp[-R 2 (k 2 + l 2 + ukl)]k a l h u c P(k)P{l). (28) 

JO JO J-1 

For a power-law initial fluctuation P(k) = k n , we obtain a final closed formula (Matsubara 
1994, Lokas et al. 1995) 

/n + 3 n + 3 3 1\ / 8\ (n + 3 n + 3 5 1\ 
^ (n) = 3F {— — 2> V ~ { n+ 7j {— — 2> V ' (29) 

where F is the Hypergeometric function. 

In the case of the top-hat filter whose Fourier transform is given by 

3 

w(kR) = — (sinkR- kRcoskR), (30) 

[fx, IX J 

the final form of S becomes very simple as follows (Juszkiewicz et al. 1993, Bernardeau 
1994) 

5 F (n) = y-(n + 3). (31) 

This formula is not only valid for pure power-law initial fluctuations but also for general 
power spectra with effective spectral index defined at the smoothing radius R as 

d In Up . . 

n = — 32 
d\nR v ; 

Equations (29) and (31) are only the leading-order contribution and more higher-order 
effects might change them considerably. Thus it is quite important to compare these 
analytic formulas with fully nonlinear numerical simulations and clarify validity of the 
perturbative formulas. There are many works on this topic and the analytic predictions 
show surprisingly good agreement with numerical simulations, even at a ~ 1 (e.g. Baugh, 
Gaztanaga & Efstathiou 1995, Hivon et al. 1995, Juszkiewicz et al. 1995, Lokas et al. 
1995). 
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So far we have discussed skewness S with fixed smoothing methods. For the third-order 
moments (Sp), the second-order effects caused by the gravitational evolution and that 
caused by the adaptive smoothing are decoupled, as we can see from equations (10) and 
(11). Thus we can write the skewness parameter for adaptive methods in the following 
forms 

S G = S F + AS G , (33) 
S s = S F + AS S . (34) 

Here AS G and ASs are the correction terms caused by the adaptive smoothing methods. 
In the next two subsections we calculated these terms explicitly. 



3.1. Gather Approach 

First we calculate the correction term AS G for the gather approach. With equation (10) 
this term is easily transformed to the following equation (see Appendix A.l for derivation) 

For a power-law spectrum we have simple equation below 

AS G (n) = (n + 3). (36) 

In derivation of equation (35) we only use Gaussianity of the one point PDF of the linear 
smoothed field 8 R . Therefore, these formulas do not depend on the shape of the smoothing 
filter nor the cosmological parameters VL or A. Furthermore they are valid also in the 
redshift space, if we use the the distant observer approximation. Thus equation (35) has 
strong predictability. 

Hivon et al. (1995) perturbatively examined the skewness parameters in redshift 
space and evaluate them both for the top-hat filter and the Gaussian filter. They also 
compared their analytic results with numerical results. They found that these two show 
agreement only in the range a < 0.1, in contrast to the skewness parameter in the real 
space a < 1.0. They commented that this limitation is mainly due to the finger of god 
effects {e.g. Davis & Peebles 1983). 

Here we use their analytic results and combine our new formula with them. In figure 
1 we present the skewness parameters for various spectral indexes n both in the real and 
redshift spaces. For simplicities we limit our analysis for the Einstein de-Sitter background. 
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For the Gaussian filter, the skewness parameter by the gather method is a increasing 
function of spectral index n both in real and redshift spaces. This dependence is contrast 
to the skewness with the fixed smoothing method. Comparing the skewness in the real and 
the redshift spaces, n dependence of the gather method is somewhat weaker in real space 
than in redshift space, but this tendency is also different from the fixed smoothing. 

For the top-hat filter, there is no spectral index n- dependence in the real space. We 
have S = 34/7 which is the same as the unsmoothed value (Peebles 1980). Bernardeau 
(1994) pointed out that the skewness S filtered with the top-hat filter in Lagrangian space 
does not depend on the power spectrum and is given by S = 34/7. As the adaptive 
smoothing is basically Lagrangian description, this fact seems reasonable. In the case of the 
redshift space we have a fitting formula below 

35.2 

Sc( n ) — — 0.15(n + 3), (Einstein de Sitter background) (37) 

which is based on formula (49) of Hivon et al (1995). Again n dependence is very weak 
and becomes weaker for < 1 (see Fig.4 of Hivon et al. 1995). Finally, we comment the 
possibility that our perturbative treatment of the redshift space skewness becomes worse 
in the adaptive methods than in the fixed smoothing method. In the adaptive methods, 
smoothing radius of a high density region becomes smaller and the (strongly nonlinear) 
finger of god effects might not be suppressed well. 



3.2. Scatter Approach 



Next we calculate the correction term AS$ for the scatter approach. We only discuss 
the real space density field smoothed with the Gaussian filter (eq.[3]). From equation (12) 
we obtain the following equation (see Appendix A.l), 



AS SR (n)a R = 2 



dk dl 

(2tt) 3 (2tt) 3 



cxp 



(3Z 2 + 2k 2 + 2k ■ l)R 2 



P(k)P(l)(k + iyR 2 . (38) 



The six dimensional integral dkdl is simplified to a three dimensional integral dkdldu as in 
equations (25) and (28). Then we have the following relation 

(3Z 2 + 2k 2 + 2klu)R 2 ' 



AS S R(n)a R 



1 



dk 



roo rl 

L *L 



duexp 



4tt 4 Jo 

xk 2 l 2 P(k)P(l)(k 2 + l 2 + 2klu)R 2 . 
For a pure-power law fluctuation we obtain the following analytic formula 



(39) 



AS s {n) = -2( n+3 )/ 2 3-(" +5 )/ 2 (n + 3)|^(n + 3)F(^ 



5 + n 5 + n 5 1 
2 ' 2 '2' 6 
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In contrast to the previous gather approach, this result is valid only to the real space 
skewness with the Gaussian filter. In table 1 we present numerical values of ASs{n). In 
figure 2 we show S,s{n) as a function of the spectral index n. We can see that n dependence 
is similar to the gather approach but now it becomes weaker. If we change n from —3 to 
— 1, skewness S changes ~ 25% for the scatter approach, ~ 45% for the gather approach, 
and ~ 38% for the fixed smoothing. 



4. Statistics of Isodensity Contour 

The genus number is a topological quantity and defined by the number of the homotopy 
classes of closed curves that may be drawn on a surface without cutting them into two 
pieces. This definition seems highly mathematical, but there are more intuitive methods to 
count the genus number. First one is to notice the number of holes and isolated regions of 
the surface in interest. Second one is to count stationary points of the surface along one 
spatial direction (Adler 1981, Bardeen et al. 1986). With these equivalent methods, we can 
calculate the genus density as follows 

„ , iV(holes) — ^(isolated regions) , 

Genus density = — ^ - (41) 

volume 

iV(maxima) + A^(minima) — iV(saddle points) 
2 x volume 

For example, in the case of one-sphere, we have iV(holes) = and ^(isolated regions) = 1 
and genus number becomes —1. We obtain the same result with equation (42). The genus 
number density of isodensity contour of the large-scale structure is a powerful measure 
to quantify connectivity of galaxy clustering, such as, filamentary networks, sheet-like or 
bubble-like structures. The genus density of a high density contours is expected to be 
negative as the surfaces would show disconnected meatball-like structure. But the genus 
density for contours around the mean density 5 ~ would be positive as they would look 
like highly connected sponge-like structure (Gott, Melott & Dickinson 1986). The genus 
density as a function of the matter density threshold is called the genus statistics and has 
been widely investigated both numerically and observationally (Gott, Weinberg & Melott 
1987, Weinberg, Gott & Melott 1987, Melott, Weinberg & Gott 1988, Gott et al. 1989, 
Park & Gott 1991, Park, Gott & da Costa 1992, Weinberg & Cole 1992, Moore et al. 1992, 
Vogeley, Park, Geller, Huchra & Gott 1994, Rhoads, Gott & Postman 1994, Matsubara & 
Suto 1996, Coles, Davies & Pearson 1996, Sahni et al. 1997, Protogeros & Weinberg 1997, 



-13- 



Coles, Pearson, Borgani, Plionis & Moscardini 1998, Canavezes et al. 1998, Springel et al. 
1998) 

Usually we use the local expression (42) to analytically study the genus statistics. For 
the genus density of isodensity contour at v = 5 /a, this expression is written as follows 
(Doroshkevich 1970, Adler 1981, Bardeen et al. 1986, Hamilton, Gott & Weinberg 1986) 

1 / ~ ~ « ~ ~ ^ 



G ( u ) = {$Drc[5(x) - iya]5 Drc [d 1 5(x)]5 Drc [d 2 S(x)]\d 3 S(x)\(d n S(x)d 22 5(x) - d 12 S(x) 

(43) 

where 5dtc{.-) represents the Dirac's delta function. The first one 5dtc[&{x) — vcr] specifies the 
contour v = 8/(7. The second and third ones 5r> rc [di5(x)}, 5r, rc [d 2 5(x)] specify the stationary 
points along x 3 direction. The term (d n S(x)d 22 S(x) — di 2 5(x) 2 ) is the determinant of the 
Hesse-matrix and assigns proper signatures for the stationary points corresponding to signs 
of equation (42). Even though equation (43) introduces a specific spatial direction (x 3 -axis), 
Seto et al. (1997) derived a rotationarilly symmetric formula, and studied nonlinear 
evolution of the genus statistics using the Zeldovich approximation (Zeldovich 1970) 

In the case of an isotropic random Gaussian fluctuation which is usually assumed as 
the initial condition of the structure formation, the complicated formula (43) is simplified 
to (Doroshkevich 1970, Adler 1981, Bardeen et al. 1986, Hamilton, Gott & Weinberg 1986) 




^) = — (^) e"" l\l - ,/). (11) 

where o\ is defined as 

a\ = ((V5) 2 ) = / dkP(k)k 4 w(kR) 2 . (45) 
x 1 2n 2 Jo 

Nonlinear evolution of the genus statistics had been studied using N-body simulations, but 
most analytical predictions for the genus statistics have been based on the linear formula 
(44). To compare it with observed distribution of galaxies we have to use sufficiently large 
smoothing radius to reduce nonlinearities. However, such a large smoothing radius is not 
statistically preferable for the finiteness of our survey volume. 

Matsubara (1994) improved this difficulty by taking into account of weakly nonlinear 
effects in the genus statistics (see also Hamilton 1988, Okun 1990, Matsubara & Yokoyama 
1996, Seto et al. 1997). He used the multidimensional Edgeworth expansion method and 
added the first-order nonlinear correction to the linear formula. His result is written as 

3/2 



a l \ ^ 2 /2 



(2tt) 2 \3a 2 



H 2 (v) + a (^H 5 (u) + —H 3 (u) + WH^v)) + 0(a 2 ) 



(46) 
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This formula is valid for statistically isotropic and homogeneous weakly random Gaussian 
fields. Here functions H n (u) = ( — l) n e v2 l 2 (d/ dv) n e~ v2 1 2 are the Hermite polynomials, 



H (x) = 1, H 1 (x)=x, H 2 (x)=x 2 -1, 

H 3 (x) = x 3 -3x, H 5 (x) = x 5 - 10x 3 + 15x. (47) 

In equation (46), S is the skewness parameter discussed in the previous section. T and U 
are called the generalized skewness parameters and defined by 

U = ~ <V<5 ■ V6A6) . (48) 

Matsubara & Suto (1996) examined the perturbative formula (46) using N-body simulations. 
For power- law spectra with n — 1, and —1, they found that this formula are in reasonable 
agreement with numerical results in the range —0.2 < vo < 0.4. 

Next let us briefly summarize the area statistics N 3 (u) which were proposed by Ryden 
(1988) and investigated detailedly by Ryden et al. (1989). The area statistics are defined as 
the mean area of isodensity contour surface per unit volume. For statistically homogeneous 
and isotropic fluctuations, the area statistics are equal to twice the mean number of 
isodensity contour crossings along a straight line of unit length. These two statistics are 
thus equivalent (beside factor 2), but the contour crossing statistics are easier to compute 
numerically (Ryden 1988). 

As in equation (43), the area statistics for isodensity contour 5 = vo is written as 

N 3 {v) = (S Drc [S(x) - ua]\V5(x)\) . 

In the case of isotropic Random Gaussian fluctuations, we have the following formula 
(Ryden 1988) 

w = H^f ^ (49) 

Weakly nonlinear effects on the area statistics can be discussed with a similar technique 
used to derive equation (46) (Matsubara 1995). In this case, we need information of the 
density field up to its first spatial derivative and nonlinear correction is expressed in terms 
of two parameters S and T as follows 



N 3 (is) 




l + ^[^H 3 ^) + ^H 1 (u))+0(a 2 ) 



(50) 



- 15 - 



In the rest of this section we consider nonlinear effects of the adaptive smoothing 
methods on the genus and the area statistics. We calculate the generalized skewness 
parameters T and U both for the gather approach and the scatter approach. We limit our 
analysis to the real space density field smoothed by Gaussian filters. 



4.1. Reparameterization of Isodensity Contour 

Equation (46) is weakly non-Gaussian genus density for isodensity contour surfaces 
parameterized by the simple definition v = 8 /a. There is another conventional method to 
name contour surfaces. In this method, we notice the volume fraction / above the density 
threshold of the contour in interest (e.g. Gott, Melott & Dickinson 1986, Gott et al. 1989), 
and parameterize the contour using value v r defined by 

u r = eri-\f), (51) 

where the suffix r indicates " reparameterization" and erf (a;) is the error function defined 
by 

erf(x) = --L J™ dye- y2/2 . (52) 

This procedure is a kind of Gaussianization. Two methods coincide v = v r when the one 
point PDF P(v) is Gaussian distributed. If we use this new parameterization, the genus 
curve is apparently invariant under a monotonic mapping of the density contrast field 
S. Furthermore, it has been long known that the genus curve with v r parameterization 
(51) nearly keeps its original symmetric shape (eq.[44]) in the course of weakly nonlinear 
gravitational evolution of density field (e.g. Springel et al. 1998 and references therein). 
Almost the same kind of tendency has been confirmed for the area statistics (Ryden et al. 
1989). Weakly nonlinear area density with u r parameterization remains at its linear shape 
very well. Here, let us relate these two parameterization methods for weakly nonlinear 
regime. Using the Edgeworth expansion method, the one point PDF P(v) is written in 
terms of the skewness S up to the first-order nonlinear correction 



2tt 



(53) 



Juszkiewicz et al. (1995) examined this approximation using N-body simulations. They 
found that the above formula is accurate until aS reaches 1. Therefore the inequality 



< 5-^0.2, (54) 
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would be a standard for the validity of the perturbative analysis in this subsection. The 
volume fraction f(u) above the threshold 5 — va is given by 



Q 

dxP(x) = erf(z/) + {e~ v \u 2 - 1)} + 0{a 2 ). 



6V27T 

With equations (51) and (55) we obtain correspondence between v and v r as follows 

V = Vr + ^{>t-l} + 0(0>). 

Finally the genus density G r (v r ) in this new parameterization is given by G(u) = G r (u, 
and written as 



(55) 



Griff, 



1 



2 \ 3/2 



(2tt) 2 \3a 2 J 
+H 1 (v)(-S + 3U))+0(a 2 ) 



H 2 (u r ) + a[H 3 (u)[-S+-T 



Similarly we have the following result for the area statistics 



N 3r (u r .) 



9 / ^ \ 1/2 

I 1 I -"' 2 



vll2 



7T \ 3(T 2 



l + a 



(56) 



(57) 



In this case, the first nonlinear correction is simply proportional to v r (see eq. [47]) and it 
completely vanishes when we have S = 3/2T. 

Later in §4.3 and §4.4, we will confirm that the nonlinear correction (proportional to 
a) for the genus and area statistics with the fixed smoothing method are very small for v r 
parameterization, as experientally known in N-body simulations. In the followings, we use 
these two parameterizations v and v r . 



4.2. Generalized Skewness 

The generalized skewness T and U are basic ingredients to perturbatively evaluate the 
weakly non-Gaussian effects on the genus and the area statistics. For the Gaussian filter 
with a fixed smoothing radius, explicit formulas valid for the power-law initial fluctuations 
were derived by Matsubara (1994). They are given as follows 



FR 



'n + 3 n + 5 3 1\ / 18\ /n + 3 n + 5 5 1 
3F — — , — — , -, - - in + — )F 



2 '2'4/ V 7/ V 2 ' 2 '2'4 
t ^m^U), (58) 



105 V 2 ' 2 '2' 4 
n /n + 5 n + 5 5 1\ 7n + 16 (n + 5 n + 5 7 1\ 

^ = ^ ^ 2' ij F l^~'^~'2'4j- (59) 
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As shown in the case of the skewness parameter S analyzed in §3.1 and §3.2, the second-order 
(first nonlinear) effects caused by the adaptive smoothing methods are decoupled from 
those induced by gravitational mode couplings. Thus we can calculate them separately and 
express the total values in forms similar to equations (33) and (34). 

First, we analyze correction terms AT GR and AU GR for the gather approach. We define 
these terms by the following equations 



T 



GR 



T FR + AT 



GR, 



u, 



GR 



u. 



FR 



AU, 



GR- 



(60) 



After some tedious algebra using equation (10), we obtain the leading-order correction 
terms as follows (Appendix A. 2) 



AT 



GR 



(n + 4), 



AU GR = g(n + 5). 



(61) 
(62) 



These results are similar to the correction term for the skewness AS GR = n + 3 (eq.[36]) 
which does not depend on the shape of the filter function. However, situation is not so 
much simple here. For the Gaussian filter we have the following relation, 

,dS R (x) 



R- 



dR 



= R 2 AJ R (x). 



(63) 



This relation plays important roles to derive equations (61) and (62). But it does not 
hold for general filter functions and the simple results given in equations (61) and (62) are 
specific to the Gaussian filter. We summarize numerical data for the parametersT^ and 
Ugr in Table 3. 

Generalized skewness for the gather approach is much more complicated. If we write 
down them in the form 



Tsr — T SR + AT SR , U SR — U S r + AU SR , 



(64) 



correction terms AT SR and AUs R are written in the manner similar to equation (39) as 
follows (Appendix A. 3) 



I rOO roo r\ 

2AT SR (n)a R a 2 R = — dk dl J ^duexp 



(3/ 2 + 2k 2 + 2klu)R 2 



-AU SR (n)a 4 1R = 



xk 2 l 2 P(k)P(l)(k 2 + l 2 + 2klu)(k 2 + l 2 + klu)R 2 , 

(3/ 2 + 2k 2 + 2klu)R 2 



] rco poo rl 

— - / dk I dl I duexp 

67T 4 Jo JO J-l 



xkH 4 P(k)P(l)(k 2 + l 2 + klu){\ - u 2 )R 2 



(65) 



(66) 
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We can calculate them explicitly as follows 



AT. 



SB 



SB 



2(n+5)/2^-{n+9)/2 



4(n + 3)F ( 



n + 5 n + 5 5 1 



V 2 ' 2 '2'6 

5. .. ,„/n + 5n + 7 5 1 
+ -(n + 3)(n + 5)F 



V 2 



2 ' 2' 6 



13(n + 5) n /n + 3 n + 7 3 1 



12(n + 



2 V 2 ' 2 ' 2' 6 

n + 5 n + 5 3 1 

V 2 ' 2 '2' 6 



3)f( 



_2("-+5)/2g-(n+5)/2 



4F 



n + 5 n + 5 5 1 



2 ' 2' 6 



5(n + 5) /n + 5 n + 7 5 1 
f 9 V 2 ' 2 '2' 6 



-4F 



n + 5 n + 5 3 1 



2 ' 2' 6 



(67) 



(68) 



We present numerical data of the parameters T$r and Usr in Table 4. Note that the 
magnitude of generalized skewness T and U becomes very small in the scatter approach. 
This fact becomes important in the next subsection. 



4.3. Weakly Nonlinear Genus Statistics 

In figures 3 to 5, we show the weakly nonlinear genus density smoothed by three 
different methods (fixed, gather and scatter), using two types of parameterizations v and 
v r . All of these curves are smoothed by the Gaussian filter (eq. [3]). We plot the normalized 
genus curves [] G(v)/G(0) or G r (v r )/G r (0) to see deviation from the symmetric linear 
genus curve oc (1 — u 2 ) exp(— z/ 2 /2). Upper panel of figure 3 is essentially same as Fig.l of 
Matsubara (1994). 

First, comparing upper and bottom panels of Fig.3, we can confirm the fact 
experientally known in N-body simulations. Weakly nonlinear genus curves for the 
fixed smoothing method are very close to the linear symmetric shape in the case of v r 
parameterization (e.g. Springel et al. 1998). For a spectral index with n ~ —1, three curves 
for a = 0, 0.2 and 0.4 are nearly degenerated. 

In Fig. 4 we present genus curves with the gather smoothing. From Figs. 3 and 4 it 
is apparent that deviations from the linear curves become larger in the gather approach, 



1 Note that the amplitude of G(0) or G r (0) are not changed by the first-order correction of a. 



-19- 



especially in tail parts. But these deviations become smaller with using parameter v r . 

Bottom panel of Fig. 4 is calculated under conditions (gather approach and parameter 
v r ) similar to Fig. 7 of Springel et al. (1998) which is obtained from N-body simulations. 
However overall shapes of these two are different. The minimum value of G r (u r ) are attained 
around the point v r ~ 1.5 in our result, but this point is v r ~ —1.5 in theirs. This difference 
might be caused by the difference of adopted filter function. We use the Gaussian filter but 
a different kernel (a spline kernel that is often used in SPH simulations) is adopted in their 
calculation (Monaghan & Lattanzio 1985). 

In figure 5 we show results for the scatter approach. Nonlinear effects are more 
prominent than two cases analyzed earlier. As shown in equations (46) and (56), nonlinear 
correction of the genus curves are written by combination of terms proportional to 
parameters S, T and U. Some of their contribution cancel out, as realized in the case of 
the fixed or gather smoothing methods. However, amplitude of parameters T and U for the 
scatter approach becomes very small (see Table 4), and cancellation becomes weaker. 

Based on the definition of the genus density (eq.[41]), nonlinear evolution of isodensity 
contours is sometimes described with such terminologies as, sponge-like (connected 
topology) or meatball like (disconnected topology). In the case of Random Gaussian initial 
fluctuations, linear theory predicts symmetry of the genus statistics with respect to the sign 
of density contrast 5, and geometry of both high and low density tails look meat-ball like 
with negative genus density. If we use the fixed or gather smoothing methods, nonlinear 
effects make the genus number of a high density contour (e.g. v — 2) smaller, and topology 
of that region becomes more meatball-like (see figures 3 and 4). In contrast, contour of low 
density threshold (e.g. v = —2) is transformed in the direction of sponge-like topology, as 
quantified by the increase of genus number density (see figures 3 and 4). 

It is not easy to understand behaviors of Fig. 5 for the scatter approach. But, changing 
point of views, we can discuss characters of this approach by comparing figures for various 
smoothing methods. To characterize nonlinear effects accompanied with the scatter 
approach, we apply the topological interpretation mentioned in the last paragraph. As 
shown in Fig. 6 where various smoothing methods are compared, the scatter approach makes 
low density regions more meatball-like. But high density regions are transformed in the 
opposite direction. This trend shows remarkable contrast to the nonlinear gravitational 
effects traced by the simple fixed smoothing method. 
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4.4. Weakly Nonlinear Area Statistics 

In Figs. 7 to 9, we plot the weakly nonlinear area statistics, using equation (50) for 
v parameterization and equation (57) for v r parameterization. As in the analysis of the 
genus statistics, we normalize amplitude of the area density by AT 3 (0) or N 3r (0). Comparing 
upper and bottom panels of figure 7, it is apparent that v r parameterization is very 
effective to keep the original linear shape against weakly nonlinear effects. This fact has 
been confirmed experimentally in N-body simulations (Ryden et al. 1989) and is quite 
similar to the situation in the genus statistics explained in previous subsection. With v r 
parameterization (bottom panel of Fig. 7), three curves for a = 0,0.2 and 0.4 are almost 
completely overlapped for all spectral indexes n. This fact seems reasonable as we have 
Sf — 3/2Tp for the fixed smoothing method (see eq.[57] and Tables 1, 3). 

If we use the adaptive smoothing methods, weakly nonlinear correction on N^u) [y 
parameterization) are considerable as shown in upper panels of Figs. 8 and 9. This correction 
becomes apparently smaller for the gather approach, but not for the scatter approach. We 
have already commented that the first-order nonlinear correction for the area statistics is 
characterized by two parameters S and T. For the scatter approach, T parameter is very 
small for spectral indexes n > — 1, and cancellation mentioned in §4.3 is not effective. 

5. Summary 

Observational analysis of galaxy clustering is one of the central issues in modern 
cosmology. Various methods have been proposed to quantify the clustering, and many of 
them (e.g. topological analyses of isodensity contour) are based on continuous smoothed 
density field. However what we can observe directly is distribution of point-like galaxies. 
Thus smoothing operation is crucially important in the analyses of the large-scale structure. 
From theoretical point of views, filters with spatially constant smoothing radii are natural 
choice and have been widely adopted so far. But we should notice that there are no strong 
convincing reasons to stick to this traditional method. 

There are few galaxies in void regions even at semi-nonlinear scales. In these regions 
density field obtained with fixed smoothing radius might be considerably affected by the 
discreteness of (point-like) mass elements, and might hamper our analyses of the cosmic 
structures. Adaptive smoothing is basically Lagrangian description, and we use nearly same 
number of particles to construct smoothed density field at each point. Thus it is quite 
possible that the adaptive methods are more efficient than the fixed methods to analyze 
the large-scale structure. Actually, Springel et al. (1998) have recently pointed out that 
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using adaptive niters, signal to noise ratio of the genus statistics is improved even at weakly 
nonlinear scales > lO/i^Mpc. 

In this article, we have developed a perturbative analysis of adaptive smoothing 
methods that are applied to quantify the large-scale structure. Even though adaptive 
methods might be promising approaches in observational cosmology, this kind of analytic 
investigation has not been done so far. Our targets are weakly nonlinear effects induced by 
two typical adaptive approaches, the gather approach and the scatter approach (Hernquist 
& Katz 1989). The concept of these methods is easily understood with equations (8) and 
(9). The gather approach is easier to handle analytically. Numerical costs dealing with 
discrete particles' systems are also lower in this approach (Springel et al. 1998). Effects 
caused by these two adaptive methods start from second-order of 5 in perturbative sense. 
They modify quantities which characterize the nonlinear mode couplings induced by gravity 
(e.g. Peebles 1980). 

In §3 we have investigated the skewness parameter S which is a fundamental measure 
to quantify asymmetry of one point PDF. We have shown that the skewness for a gather 
top-hat filter does not depend on the spectral index n in real space, and very weakly 
depends on it (S ~ 35.2 — 0.15ra: Einstein de-Sitter background) in redshift space. In 
the case of Gaussian filter, the skewness parameters show similar behaviors both in the 
scatter and gather approaches. They are increasing functions of n, in contrast to the fixed 
smoothing method. 

Next in §4, the genus and area statistics have been studied with Gaussian adaptive 
filters. Our analysis is based on the multidimensional Edgeworth expansion explored by 
Matsubara (1994). We use two quantities v{= S/a) and v r to parameterize isodensity 
contours. The latter v r is defined by the volume fraction above a given density threshold 
(Gott, Melott, & Dickinson 1986). It is explicitly shown that using this parameterization, 
two statistics with the fixed smoothing method are very weakly affected by semi-nonlinear 
gravitational dynamics, as experientally confirmed by N-body simulations. For the gather 
smoothing, we found that the v r - parameterization is more effective to keep original linear 
shape of the area statistics than of the genus statistics. 

The parameters S, T and U which characterize the nonlinear corrections of isodensity 
contour depend largely on the filtering methods. We can characterize nonlinear effects of 
these methods in somewhat intuitive manner, using results for the genus statistics. The 
scatter approach makes low density tails more meatball-like, but high density tails are 
transformed in the direction of sponge-like (connected) topology. This is a remarkable 
difference from fixed or gather smoothing methods. 
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Our investigation in this article has been fully analytical one, using perturbative 
technique of cosmological density field. Numerical analyses based on N-body simulations 
would play complementary roles to results obtained here, and thus are very important. 
Apart from numerical investigations, perturbative treatment given in this article would be 
also developed further in several ways. The smoothed velocity field is crucially important 
material in observational cosmology as it is supposed to be less contaminated by effects of 
biasing (Dekel 1994, Strauss & Willick 1995). But our observational information is limited 
to the line of sight peculiar velocities only at points where astrophysical objects exist. 
Thus construction of smoothed velocity field contains similar characters as discussed in 
this article (e.g. Bernardeau & van de Weygaert 1996). There is another (more technical) 
problem that has not mentioned so far. In this article we have only studied spherically 
symmetric filter functions. Springel et al. (1998) have shown that signal to noise ratio of 
the genus curves is further improved by using a triaxial kernel, taking account of tensor 
information of local density field. This point must be also worth studying. 

The author would like to thank J. Yokoyama for discussion and an anonymous referee 
for useful comments. He also thanks H. Sato and N. Sugiyama for their continuous 
encouragement. This work was partially supported by the Japanese Grant in Aid for 
Science Research Fund of the Ministry of Education, Science, Sports and Culture No. 3161. 
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A. Derivations of Parameters 

In this appendix we derive expressions for the correction terms AS, AT and AU given 
in the main text. First we perturbatively expand the density contrast field smoothed by an 
adaptive filter as 

5 A (x) = 5i(a?) + 5 2 (x) + 5 2A (x) + ■■■, (Al) 

where Si(x) is the linear mode, 8 2 (x) is the second-order mode induced by gravity, and 
5 2A is the second-order correction term caused by an adaptive smoothing (the suffix A 
represents "adaptive"). Then the third-order moment for S A (x) is given as 



[6%) = 3 ( 6 2 6i) + 3 {S 2A Si) + ■■■. (A2) 

Thus the first-order correction term for the third-order moment is given as 

3(yj). (A3) 

In the same manner we have the following correction terms for (5 A V 2 5 A ) and 
(V5 A ■ V5 a V 2 <5 A ) as 

[S 2 1 V 2 5 2A ) + 2(5 1 5 2A V 2 5 1 ), (A4) 
, V5 1 ■ V5 1 V 2 5 2A ) + 2 (V5iV5 2 aV 2 5i) . (A5) 



We can write down the second-order correction terms 5 2A for the gather and scatter 
approaches with smoothing radius R (see eqs. [10] and [11]) 

1 d 1 d 

$2Gr{x) = --5 1R {x)R—5 1R {x) + -R—a 2 R , (A6) 

<W*) = -fjdx'd R W(\x'-x\;R)5 1 (x')5 1R (x'). (A7) 

where S±r(x) represents the smoothed linear mode, W^dcc' — x\; R) is a filter function. The 
variance a R of the matter fluctuations is given as 

a 2 R = (5 2 R ) + 0(5*) (A8) 

= j ^Y^ kR ) 2p ^) + ^ ( A9 ) 

where w(kR) is a Fourier transformed filter function. 

For the Gaussian filter (see eq. [3]) the above equations are given with the linear Fourier 
modes S 1 (k) as 



. , f dk dl 

hGR{x) = J ( 2? r)3 ( 27 r)3 6XP 



(I 2 + k 2 )R 2 



k 2 R 2 

«Ji(fc)«Ji(l)— — exp[i(fe + I) ■ x] 
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1 d 9 

. . . f dk dl 

62SR(X) = -Jj^fJ^f^ 
x exp[i(k + I) ■ x], 

Next we comment on the ensemble average of variables. We assume that the linear 
Fourier modes of density fluctuation are random Gaussian distributed. If variables 
{A, B, C, D} obeys multivariable Gaussian distribution, we generally have the following 
relation 

(ABCD) = (AB) (CD) + (AC) (BD) + (AD) (BC) . (A12) 
For the linear Fourier modes the above equation becomes 

(<5i(fc)<Si(0<Si(™)<*i(«)> = (27i) 6 P(k)P(l)5 Drc (k + m)5 Drc (l + n) 

+ (2nfP(k)P(m)5 Drc (k + l)5 Drc (m + n) 

+ (2ixfP(k)P(l)5 Drc (k + n)5 Drc (l + m). (A13) 

Here 5dtc( ) is the Dirac's delta function and P(k) is the matter power spectrum. We 
evaluate expressions (A3)-(A5) using relations (A12)-(A13). 



(A10) 



(2Z 2 + k 2 + 2k ■ l)R 2 



<$i(fc)*i(0 



(k + lfR 2 



(AH) 



A.l. Skewness 

For the gather approach the real-space representation (A6) is more convenient. Using 
property (A12) we obtain the following result 

= ( SiR m Si, <) + r 2 «jR a « (A15) 

= (Ai6) 

The correction term for the skewness Sq is written as equation (35) 

A5g = = "^d/T* = - dhif ■ (A17) 

For power-law models we have a simple relation a\ oc R~ n ~ 3 , and the above expression 
becomes 

A,S G = (n + 3). (A18) 
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Here we should notice that relations (A17) and (A18) do not depend on the choice of filter 
functions. 

Next let us evaluate the correction term for the skewness parameter in the case of the 
scatter approach. In this case we limit our analysis for a Gaussian filter (eq. [3]). From 
equation (All) we have 

i , . „ . . 9 \ r dk dl dm dn .. 9 9 

(2l 2 + k 2 + m 2 + n 2 + 2k ■ l)R 2 

exp 



x (5i(k)5i(l)5i(m)5i(n)) exp[i(k + I + m + n) ■ x] 



Using equation (A13) we can simplify the above integral as 



dk dl 

(2vr) 3 (2tt) 3 



exp 



(3l 2 + 2k 2 + 2k ■ l)R 2 



P(k)P(l)(k + l) 2 R 2 



(A19) 



(A20) 



Note that the integrad of this expression depends only on the information of the shape of 
the triangle determined by two vectors k and I. This triangle is characterized by three 
quantities, namely, two sides k = \k\, I = \l\ and cosine between them u = k ■ l/(kl) with 
— 1 < u < 1. We change variables from {k, 1} to {k,l,u}. The volume element is deformed 

as 

dkdl 8n 2 dkdldu. (A21) 

Thus we obtain equation (39). This equation looks somewhat complicated. For power-law 
models, however, we can easily evaluate it using mathematica (Wolfram 1996) and finally 
obtain analytical expression (40) given in the main text. 



A. 2. Generalized Skewness for the Gather Approach 

For this approach we have the following relation for a Gaussian filter 

R-^5 1R (x)=R 2 V 2 5 1R (x). 



(A22) 



Therefore the correction terms (A4) and (A5) can be written by combinations of the 
following five variables 

{5, V5, V 2 5, V 3 5, V 4 <5}. (A23) 



2 In this subsection we denote 6ir(x) simply by S. 
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For example, equation (A4) is written as 
1 



"I* 



5 2 V 2 (5V 2 5)) + 2 (5 2 (V 2 5)(V 2 5) ) - 2 (5V 2 5) (5V 2 5 



Using property (A12), the above expression is deformed as 
1 



-R 1 



3 ( V 2 5V 2 5) (55) + 2 ( V5V"5) (55) + 3 (5V 4 5) (55) + 4 (5V 2 5 



The moments appeared in the above equation can be written in terms of P(k) as 

dk 



- (5V 2 5 
V 2 5V 2 5 



(W5W5) = 
-(V 3 5V5 



2tt 2 



k*P(k)e 



V 4 M 



dk 



2tt< 



k 6 P(k)e 



-k 2 H 2 



For a power-law models (P(k) oc k n ). These integrals are evaluated respectively as 



• ,. i n ~ [i >\ a 2 r -a f n + 3\ + 5 



<R 



7 U R J 



(A24) 

(A25) 

(A26) 
(A27) 

(A28) 



With the definition of T parameter (eq. [48]) we obtain the final result that is given in 
equation (61) as 

AT GR = ^(n + A) (A29) 

To calculate the correction term (A5), let us use the Fourier space representation 
(AlO).Q It is straightforward to get 



Wi ■ V5 1 V 2 5 2A ) + 2 (V5 l V5 2A V 2 5 1 



dk dl dm dn k 2 R 2 

(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 3 

" (I 2 + k 2 + m 2 + n 2 )R 2 
xexp 

x (5 1 (k)5 1 (l)5 1 (m)5 1 (n)) 
x exp[z(fc + I + m + n) ■ x], 
x[-(m ■ n)(k + I) 2 
- (m • (fe + Z)™ 2 ) - (n • (fe + Z)m 2 ]QA30) 



3 We obtain the same result starting from equation (A6). 
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With equation (A13), the above integral becomes 

~\ I (2n)H2ny k2R2 6XP H 1 ' + ^ P ^l)[^f - 4(fe • if] 
= -/ dfc / dZ / dwexp -(/ 2 + A; 2 )i? 2 k 6 l A P(k)P(l)(l - u 2 )R 2 

67T 4 JO JO i-1 L J 

= -9^4 I k 6 P(k)exp[-k 2 R 2 ] J dll 4 P(l)exp[-l 2 R 2 ]. 

For power-law models this expression becomes (see eqs.[A26]-[A28]) 

'n + 3\ 2 /n + 5 s 



9 



4 



Using definition of C/ parameter (eq.[48]) we obtain equation (62) as 

AU GR = ^(n + 5). 



(A31) 



(A32) 



(A33) 



Note that results (A29) and (A33) are not valid for general filters. Equation (A22) 
that holds for the Gaussian filter plays crucial roles to derive them. 



A. 3. Generalized Skewness for the Scatter Approach 



First we evaluate the correction term given in equation (A4). With the Fourier space 
representation (A10) we obtain the following equation 

dk dl dm dn (k + l) 2 R 2 



5 2 l V 2 5 2A ) + 2(8 1 8 2A V 2 5 l 



(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 3 

(2/ 2 + k 2 + m 2 + n 2 + 2k ■ l)R 2 



x exp 



x (5i(k)5i(l)5i(m)Si(n)) exp[i(fe + I + m + n) ■ x], 
x[(k + I) 2 + m 2 + n 2 ] (A34) 

With the formula (A13), this twelfth-dimensional integral becomes 



2 r dk dl 

3 J (2tt) 3 (2tt) 3 6XP 



(3Z 2 + 2k 2 + 2k ■ l)R 2 



P(k)P(l)(k + lYR\k+iy+l 2 + k 2 }. (A35) 



changing variables from dkdl to dkdldu as shown in relation (A21), we obtain the result 
essentially same as equation (65) as 

(3/ 2 + 2A; 2 + 2klu)R 2 

6tt 4 



. . . I roo poo rl 

6 2 V 2 5 2A ) + 2(5 1 5 2A V 2 5 1 ) = — J dk J dl J ^duexp 



xk 2 l 2 P{k)P{l){k 2 + l 2 + 2klu){k 2 + I 2 + klu)R 2 , (A36) 
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As in the case of skewness parameter, we can evaluate this complicated integrals with 
mathematica and obtain equation (67). 

In the same manner, equation (A5) is written as 

„ r \ „/ orr7r f dk dl dm dn (k + l) 2 R 2 



(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 3 

(2Z 2 + k 2 + m 2 + n 2 + 2k ■ l)R 2 
xexp 

x (5 1 (k)5 1 (l)5 1 (m)5 1 (n)) 
x exp[i(k + I + m + n) • x], 
x[-(m • n)(k + I) 2 
- (m ■ (k + Z)n 2 ) - (n ■ (k + l)m 2 )].(A37) 



This expression is simplified to the following form 



dk dl , m \ n ^ o 
; — — — (k + lYR 2 exp 
3 J (2tt) 3 (2tt) 3V ; P 

xP(A;)P(/)[4A; 2 / 2 -4(fc-Z) 2 ]. 



(2i 2 + k 2 + m 2 + n 2 + 2k- l)R 2 



(A38) 



Again, changing variables, we obtain the expression (66) as 

1 



V5i • V5i V Z 5 2A ) + 2 ( V5 l V5 2A V z 5 l ) = 



6vr 4 
x exp 



1 - u 2 )R 2 



roc rOO pi 

j dk j dl J ^ du{ 

(3/ 2 + 2A; 2 + 2klu)R 2 



2 

.,2 , ;2 



xkH 4 P(k)P(l)(k 2 + l 2 + 2klu). (A39) 

For power-law models we can evaluate this integral using mathematica and obtain equation 
(68). 
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TABLE 1 

SKEWNESS FOR THE GATHER APPROACH (GAUSSIAN FILTER) 



spectral index n 


1 





-1 


-2 


o 
-o 


S F (n) 


3.029 


3.144 


3.468 


4.022 


4.O0 1 


AS G (n) 


4.000 


3.000 


2.000 


1.000 


n 
u 


S G (n) 


7.029 


6.144 


5.468 


5.022 


4.857 






TABLE 2 








SKEWNESS FOR THE SCATTER APPROACH 




spectral index n 


1 





-1 


-2 


o 
-o 


S F (n) 


3.029 


3.144 


3.468 


4.022 


4.O0 ( 


A s S(n) 


3.031 


2.576 


2.045 


1.277 


n 
u 


Ss(n) 


6.060 


5.720 


5.513 


5.299 


4.o0 ( 






TABLE 3 








GENERALIZED SKEWNESS FOR THE 


GATHER APPROACH 




spectral index n 


1 





-1 


-2 


-3 


T F (n) 


2.020 


2.096 


2.312 


2.681 


3.238 


AT G (n) 


3.333 


2.667 


2.000 


1.333 


0.667 


T G {n) 


5.353 


4.763 


4.312 


4.014 


3.905 


U G (n) 


1.431 


1.292 


1.227 


1.222 


1.272 


AU G (n) 


2.000 


1.667 


1.333 


1.000 


0.667 


U G {n) 


3.431 


2.959 


2.560 


2.222 


1.929 






TABLE 4 








GENERALIZED SKEWNESS FOR THE SCATTER APPROACH 


spectral index n 


1 





-1 


-2 


-3 


T F {n) 


2.020 


2.096 


2.312 


2.681 


3.238 


AT s (n) 


-2.082 


-1.908 


-1.723 


-1.451 


-0.963 


T s (n) 


-0.0623 


0.1882 


0.5892 


1.230 


2.275 


U F (n) 


1.431 


1.292 


1.227 


1.222 


1.272 


AU s (n) 


-1.265 


-1.145 


-1.027 


-0.8916 


-0.7105 


U s (n) 


0.1662 


0.1474 


0.2000 


0.3301 


0.5611 
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Fig. 1. — Skewness for the gather approach in real and redshift spaces. We use two kinds of 
filters (Gaussian and top-hat filters). The dashed-lines represent results for the traditional 
fixed smoothing method and the solid lines for the gather approach. Numerical data in the 
redshift space are based on Table 2 and 3 of Hivon et al. (1995). 



Fig. 2 




Fig. 2. — Real space skewness for the scatter approach with the Gaussian filters. The 
dashed-line represents results for the fixed smoothing method and the solid line for the 
scatter approach. 
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Fig. 3. — Normalized genus density for the fixed Gaussian smoothing (see Matsubara 1994). 
Upper panel corresponds to v- parameterization and lower to v r - parameterization. The 
solid curves represent linear curves (eq.[45]). Dotted-lines, dashed-lines show a = 0.2, and 
a = 0.4 respectively. 
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Fig. 4. — Same as Fig. 3 but with the gather approach. 




Fig. 5. — Same as Fig. 3 but with the scatter approach. 
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Fig. 6. — Nonlinear Genus curves for various smoothing methods. The solid line corresponds 
to the linear analysis, dotted to the fixed smoothing, short-dased to the gather approach, 
and long-dased to the scatter approach. We fix the spectral index at n = and nonlinearity 
at a = 0.2. 
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Fig.7 




-2 2 -2 2 



Fig. 7. — Normalized area statistics for the fixed Gaussian smoothing. Upper panel 
corresponds to v parameterization and lower to v r parameterization. The solid curves 
represent the linear prediction (eq. [45]). Dotted-lines, dashed-lines show a = 0.2, and a = 0.4 
respectively. 
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Fig. 8. — Same as Fig. 6 but for the gather approach. 
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Fig. 9. — Same as Fig. 6 but for the scatter approach. 



